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Abstract 

To extend previous results on the late time behavior of massive helds, for the Dirac field propagating 
in the D-dimensional Minkowski spacetime we calculate analytically its asymptotic tails. We find that 
the massive Dirac held has an oscillatory inverse power law tail. The frequency of the oscillations depends 
on the mass of the held and the power law decay rate depends on the dimension of the spacetime and the 
mode number of the angular eigenvalues. We also compare with previous results in curved spacetimes. 
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1 Introduction 

At present time, it is well established that the evolution of a test field around a black hole has three stages 
[T], m, [3]: a) Prompt response, b) Quasinormal ringing, c) Asymptotic tail decay. Furthermore, it is well 
known that each stage is associated to a different part of the contour in the inversion formula for the Green 
function. Thus we know that the quasinormal ringing comes from the poles, the asymptotic tail decay is 
associated with the branch cut, and the prompt response comes from the semicircle with large radius DP: 0: 

El- 

Beginning with the work by Price i: 0 has been interest in determining the asymptotic tails in several 
spacetimes m, 0,0,0- The previous studies are focused on massless bosonic fields, as the gravitational, 
Klein-Gordon, and electromagnetic perturbations 0-0- It is well known that the late time behavior of the 
massive fields is different from that of the massless fields and for example the massive Klein Gordon field has 
tails in Minkowski spacetime [5] since the different frequencies forming a massive wave packet propagate with 
different phase velocities. For massive fields propagating in curved spacetimes we find that the computations 
of the late time tails are focused on the Klein-Gordon field (see 0-[Il]), although the asymptotic tails of 
the massive Dirac field have been calculated, mainly in four-dimensional spacetimes [in]-[in] (see [H] for the 
massive spin 2 field and [22] for the massive vector field). Furthermore we know calculations of the late time 
behavior for massless fields [23] and for the massive Klein-Gordon field [23] in higher dimensional spacetimes. 

Also the different aspects of fermion’s propagation in spacetimes have not been extensively explored, as 
for the bosonic fields. Thus it is convenient to study in detail the propagation of fermionic fields in different 
spacetimes. Therefore, in this work we analytically calculate the asymptotic tails of the massive Dirac field 
propagating in the D-dimensional Minkowski spacetime (D > 4), in order to extend the results that appear 
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in Refs. [l5]-[20] and in the previous studies on the propagation of the Dirac field in higher dimensional 

spacetimes m-m- 

Owing to the problem of analyzing the dynamics of massive fields is more challenging than the correspond¬ 
ing problem for massless fields we calculate the late time tails of the massive Dirac field in D-dimensional 
Minkowski spacetime since we like to know how the initial configuration of this held evolves in the time and in 
the previous results we hnd that as t —>■ oo a massive held may behave qualitatively different from a massless 
held. Also we analyze this physical situation since calculating the asymptotic time tails of the massive Dirac 
held in D-dimensional Minkowski spacetime allows us to separate the effects of held’s mass and the curvature 
of spacetime. Thus we can study only the effect of the mass on the asymptotic time behavior, that is, this 
setting allows us to clarify what effects are caused by the mass of the held in the time evolution when the 
spacetime curvature is absent. 

For massless helds we know that the power law fall-off is sensitive to the fact that the spacetime dimension 
is odd or even [23], [32]. We think that it is relevant to determine whether the mathematical form of the 
asymptotic tails of the massive Dirac held depend on whether D is odd or even. Furthermore our computation 
allows us to analyze in a simple setting the inhuence of the spacetime dimension in the late time behavior 
of this held. Therefore, the computation of these characteristics about the late time tails of the massive 
Dirac held extend our knowledge on its behavior and can be useful when more complex higher dimensional 
spacetimes will be studied, since we expect a richer phenomenology in D > 4 and our analysis can help us 
to understand what aspects of the relaxation phenomenon depend on the existence of a body that produces 
curvature and what properties depend on the spacetime dimension and the mass of the held. 

We organize the rest of this paper as follows. In Sect. [2] we calculate the late time tails of the massive 
Dirac held propagating in the D-dimensional Minkowski spacetime. First we determine the late time tails 
for D even in Subsect. 12. H and then for D odd in Subsect. 12.21 Finally in Sect.|3]we discuss our main results. 


2 Tails of the Dirac field 


In the coordinates with i = 1,2,..., D — 2, a D-dimensional maximally symmetric spacetime has a 

line element 

ds^ = E{rfdi^ - Q{rfdr^ - H{rfdE\,_^, (1) 

where E, Q, E[ are functions of the coordinate r and dE|j _2 is the line element of the (D — 2)-dimensional 
invariant base manifold with coordinates ipi. To study the behavior of the spin 1/2 held in spacetimes of the 
form 01) we use that the Dirac equation if/ip = rmp simplihes to the coupled system of two partial diherential 
equations |33]-[35] (see also Ref. m) 


E 


dttpi + -p^dr-ipi = — *^77 + imE ) 'tjj2, 


Q 


E 


E 


Q 


H 

E 


dt4’2 — -p:dr4’2 = 'iK— — imE ) V'l, 


H 


( 2 ) 


where '0i > '02 are the components of a two spinor and k are the eigenvalues of the Dirac operator on the base 
manifold with line element dE|)_ 2 . 

For the D-dimensional Minkowski spacetime, whose metric can be written in the form o with E = Q — 1, 
H = r, and (iE |,_2 being the line element of the (D — 2)-dimensional sphere, we obtain that the system of 
equations 0 transforms into 


dti’2 - dr'ip2 




5*01 + 5^01 


02 , 


(3) 


with K = ±*(n+(D—2)/2) representing the eigenvalues of the Dirac operator on the (D—2)-dimensional sphere 
where the parameter n takes the values n = 0,1, 2,..., m- In what follows we choose k = i{n+ (D — 2)/2). 
We think that similar conclusions are obtained for the other eigenvalues. 

Dehning Ei and E 2 by 


T’l = 01 -I- 02, D 2 = 01 - 02, 


(4) 


2 


from Eqs. m we find that these functions satisfy 


A 


dr - — j F 2 = -{im + dt)Fi, ( a,. + — ) = {im - dt)F 2 . 


IK 


The differential operators that appear in the previous equations satisfy 


IK 


dr ^ — {Tim — dt) ^ = (=F^^ ~ dt) dr =b 


IK 




(5) 

( 6 ) 


Taking into account this fact, from Eqs. (O we can obtain decoupled equations for the functions Fi and ^ 2 . 
The equations take the common forrro 

= (7) 

with 

L = n+— —2, {L = n+— —, (8) 

for the function Fi (F 2 ). Furthermore we point out that in both cases the quantity L is an integer for D 
even and it is a half-integer for D odd. 

To solve Eq. d?]) we use the Green function G{r,r';t) that satishes [5], [7] 

d^G - d^G - +rr/^G = 5{t - t')5{r - r'). (9) 


The causality condition imposes that for t < 0, G{r,r';t) = 0. To find Green’s function we follow the usual 
path and use its Fourier transform 


G{r,r'-u}) = / G{r,r';t)e^‘^^dt, 


( 10 ) 


that is a solution of the differential equation 

(PG 


dP 


+ ) G = 5{r — r'). 


As is well known, G{r,r'-,oj) is given by [6], [7] 


G{r,r';u}) = 


fair',^)9bir,u;) 

W{dj) 


( 11 ) 


( 12 ) 


where fa{r',u)) {gb{r,ui)) is the solution to the homogeneous form of Eq. (fTTI) that satisfies the left (right) 
boundary condition and W(u!) is the Wronskian of these solutions. To determine G(r,r';t) from G{r,r']ui) 
we use the inversion formula m, m 


G{r,r';t) = 


^+00 + 2(5 


27r 


— oo+i(5 


G(r, r';a;)e-*"’‘dw. 


(13) 


where d is a small positive constant. 

To solve the homogeneous form of Eq. (HH) we propose that its solutions take the form [9], [ID] 


F = r^+ie-“’'(/)(r) 


(14) 


^In what follows F represents to Fi or F2. 
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( 15 ) 


with Lo = y/m? — to find that the function (j) must be a solution of the differential equation 

+ 1) - z] ^ - (L + 1).^ = 0, 

where we use the variable z = 2uir. The previous equation is of confluent hypergeometric type [38], [39] 

(Ph , ^dh , , , 

"5?+ 1^-'I* 

with parameters a = L + 1, and 7 = 2 (L + 1) = 2 a. 

In H-dimensional Minkowski spacetime to calculate the tails of the massive Dirac field we impose the 
boundary conditions 

i) The field goes to zero at r = 0 

ii) The field is purely outgoing (decaying) for w > m (m > w) 

Considering the behavior in different regions of the solutions to the confluent hypergeometric differential 
equation we find that the function satisfying i) is 


/a(r,w) = Ar^+^e ‘^’'iJ^i(a, 7 ; 2wr), (17) 

where iFi(a, 7 ;z) denotes the confluent hypergeometric function [35], [31], and the function that satisfies ii) 
is equal to 

gb{r,u}) = Sr^+^e“‘^’'f/(a, 7 ; 2 wr), (18) 

where U{a,^;z) is the Tricomi solution of the confluent hypergeometric equation [31]. In (fT7)) and (fT511 the 
quantities A and B are constants. Using the expression for the Wronskian of the special functions 1 ^ 1 ( 0 , 7 ; z) 
and C/(a, 7 ;z) ((13.2.34) of Ref. [H]), we obtain that the Wronskian of the solutions fa{r,ui) and gb{r,Cd) is 
equal tcH 


W{fa{r,uj),gb{r,uj)) = W{uj) = 


AR(2L + l)r(L + 1/2) 


(19) 


As we noted previously, for D even the parameter L is an integer, whereas for D odd the parameter L 
is a half integer. As a consequence the function gb{r,u}) behaves in a different way for D even and D odd. 
Therefore in what follows we study both cases. 


2.1 D even 

For D even the quantity L is an integer, hence the parameters a and 7 are integers, thus the Tricomi function 
that appears in the expression (IT 8 )) for gb{r,ui) takes the form [31] 


U{a, 7 ; 2wr) = 


L 

E 

s =0 


(a)s 


1 


( 2 wr) 


L+l+s’ 


( 20 ) 


and the function gb{r,ui) has a branch cut between w = —m and uj = m. Therefore the function G{r,r';u}) 
has a branch cut in the interval w G (—m, m). 

As usual, to calculate the integral of the formula (IT3t we deform the contour of integration. Thus we 
close the path of integration with a semicircle of large radius in the lower complex w plane. For t > r + r' the 
integral along this semicircle goes to zero and in our case this condition is satisfied since we are interested 
in the behavior of the Dirac field as t —>■ 00 . We notice that in our problem the function G{r,r';uj) does 
not have poles in the complex plane w. Furthermore it is well known that the late time tail of the field is 
determined by the branch cut in the integral (IT31) . To calculate this integral for D even we use the contour of 

^As usual we denote the gamma function by r(z). 
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Figure 1: Path used to calculate the integral along the cut for D even. We represent the cut by a broken 
line. 


Figure [T] with the branch of uj indicated in this Figure. We obtain that the integral along the cut simplifies 
to [5]. [TU] 


G(r, r';t) 


1 

27r 



fa{r',uje^^)gb{r, uje^^) 
VF((I;e®’^) 

W{Co) 


du! 


( 21 ) 


since the integrals along the small circles Ci and C 2 of Figured] go to zero as their radii go to zero. 

For the following computations we need to know the analytical properties of fa{r,Lu), gb{r,Lu), and W{lu) 
under a rotation z —>■ ze®’^ of the argument. From m we find that the function fa{r,uj) fulfills fa{r,uje^'^) = 
fa{r,Cd), that is, this function is single valued. To evaluate gb{r,uje‘'^) we use that U{a,j-,z) is given by [40] 


U{a,-f,z) 


TT f i-Fi(a,7;z) 

sin( 7 r 7 ) ( r(l + a — 7)r(7) 


zi-T'iFi(l + a-7,2 

r(a)r(2-7) 



( 22 ) 


and therefore for D even {L integer) we obtain that the function gb{r,ui) satisfies 


gb{r,u}e"^) 


2B (-l)^+ir(L+ 1) 
~A r(2(L + l)) 


fa{r,uj) - gb{r,uj), 


(23) 


and hence the Wronskian fulfills VF(a;e®’®) = — W{uj). 
From these results we find that (l?T]l transforms into 


G(r, r'\t) 



fair',U}) 


gbjr, tue®®®) 
kF(we®®®) 


9 b{r,uj) 

WiCo) 




_Ml!_ 

A222ir(L + |)r(L + i)(2L + 1) 


fair',uj)fair,uj)u} 




(24) 


To determine the asymptotic time tail of the Dirac field we must evaluate the previous integral in the limit 
wr <C 1 i, m- To make this calculation we use the relationship between the confluent hypergeometric 
function and the modified Bessel function of first kind Iu{z) ((13.6.9) of [39]) and in contrast to the previous 
works we consider the first terms of the expansion for I,j{z) near z = 0. From these facts we obtain that in 
the limit 1/m the previous integral is equal to 




A^riL + 


tL+ 3/2 


_F(L + 3/2) 


cos(xi) + - 


1 r 4(L + 1)^ — 1 (r^ + (r')^)r 


t { 8F(L + 3/2) r(L + 3/2) 


sin(xi) 


( 25 ) 
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where 

Xi = mt- + 0 TT. (26) 

To calculate the previous integral we use the integral representation of the Bessel function ((9.1.20) of [35]) 
and the Hankel asymptotic expansions of the Bessel function ((9.2.5) of [35]). Thus in the limit t ^ 1/m the 
contribution of the cut to the Green function is equal to 


G(r, r';t) 


(-l)^r(L + 3/2)m^-i/2 (rr')^+^ 
2^-3/2r(L +1/2)(2L+1) 


m 

r(L + 3/2) 


cos(xi) 


1 f 4(L+l)^-l 

t { 8r(L + 3/2) 


(j,2 _|_ 

r(L + 3/2) 


sin(xi) + • ■ • 


(27) 


Notice that we calculate the first two terms of the expansion for G{r,r'\t) in inverse powers of t. In 
a similar way we can compute additional terms. Therefore from the previous expression, as t —>■ oo the 
dominant contribution to the function G{r,r';t) is 


G(r, r'-t) 


2i-i(2L + l)r(L+i) t^+3/2 


(28) 


Thus in this limit the function G{r,r';t) behaves as 


G(r, r';t) 


cos(xi) 

fL+3/2 • 


(29) 


2.2 D odd 


We remind that for D odd the parameter L is a half-integer and from the expressions for a and 7 we obtain 
that the quantity 7 is an integer and the mrameter a is a half-integer. Therefore the Tricomi function that 
appears in the function gt{r,u)) is equal tco [39] 


U{a, 7; 2ujr) 


_ 1 (^ ~ 1)!(~^ + ^)2L-|-l-fe 

“ r(L+l) ^ ( 2 L + l-fc)! 

(_l)2 L-^2 ^ (Q;)fc(2d;r)^ 

(2L -t l)!r(a - 2L - 1) ^ (2L -t 2 )kk\ 

-■0(1 + k) - 0(2L -I- 1 -t- fc -b 1)). 


1 

(2wr)^ 

(ln(2wr) + tp{a + k) 


(30) 


As before, for D odd the function gh{r,u}) has branch points at a; = —m, ui = m. Furthermore the 
logarithmic term implies that the function gh{r,Cd) has an additional branch point as w —> 00. Thus we can 
choose the branch cut as in Figure [2] and using the contour of this Figure we calculate the integral along the 
branch cut. We notice that the integrals along the two paths parallel to the negative imaginary axis have a 
sum equal to zero and following a similar method to that for D even we find that for calculating the integral 
along the cut we must compute an expression similar to that given in m- To calculate this integral for D 
odd we must determine the analytical properties of the functions fa{r,uj), gi,(r,Lu), and W{uj). For D odd we 
hnd that the function fa(r,uj) is single valued, as for D even, but we find that the function gb{r,uj) fulfills 

gbir,uje^^) = - -b 3/2) 


and hence W(u}) transforms in the form bF(a;e*’^) = W(u}). 

From these facts we find that for D odd the analogous of Eq. dMj) is 


G{r,r';t) = 


i(_l)^+l/22-(2L-Hl) r" 

A2r(L + |)r(L + i)(2L + i)y_, 


fair',Uj)fair,UJ)U} 




^As usual we denote the logarithmic derivative of the gamma function by 


(32) 
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Figure 2: Path used to calculate the integral along the cut for D odd. We represent the cut by a broken line. 


Using a similar method to that for D even and taking into account the integral pK]) we obtain that the 
integral along the cut produces that as t ^ oo, the function G{r,r']t) behaves as 


G{r,r';t) = 


i(—+ 3/2 )to-^ 

2L-i/2r(L + i/ 2)(2L + 1) t^+3/2 [r(L + 3/2) 


■ cos(xi 


+ - 


1 r 4(L + 1)^ — 1 (r^ + (r')^)r 


sin(xi) 


t [ 8r(L + 3/2) r(L + 3/2) 

Therefore for D odd the dominant term in the late time tail of the Dirac field is 


G{r,r';t) 


■cos(xi), 


2 L-i/ 2 ( 2 L+l)r(L +i) t^+3/2 

and hence in the limit t ^ l/m, the function G{r,r';t) behaves as in the formula 


(33) 


(34) 


3 Discussion 


in 


the form [B], [7] 


As is well known, using G{r, r']t) we write the solution to Eq. ([7]) ^ 

F(r,t) = J dr'G(r, r'\t)dtF{r',Q) + J dr'dtG{r, r';t)F{r\0) 
D even or D odd and initial data satisfvine F(r' O') = 0 or ( 


(35) 


therefore we find that for D even or D odd and initial data satisfying F{r', 0) = 0 or dtF{r', 0) = 0, the Dirac 
field propagating in the Minkowski spacetime has an oscillatory inverse power law tail of the form 


for the first case or 


F{r,t) - 

cos(xi) 

(36) 

tL+3/2 ’ 

F{r,t) - 

sm(xi) 

(37) 

tL+3/2 ’ 


or the second case. Thus in the D-dimensional Minkowski spacetime as t —>■ oo the Dirac field uecays witn 
i power t“(^+3/2) and oscillates in the form cos(xi) or sin(xi). We point out that these oscillations have a 
)eriod 2tt jm and this period is only determined by the value of the mass. 

It is convenient to notice that for initial data given by dtF{r',0) ^ 0 and F{r',0) ^ 0, for 
)ehavior of the Dirac field is 


a power and oscillates in the form cos(xi) or sin(xi). We point out 

period 2tt jm and this period is only determined by the vah 
It is convenient to notice that for initial data given by 
behavior of the Dirac field is ^ 

[gcos(xi)-f mgsin(xi)] 

with 


the Dirac field decays with 

-- ^ 

t ^ oo the 


d = J dr'{r')^^^dtF{r' ,0), J 


(38) 

(39) 
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This is due to the existence of the factor cos(xi) in the asymptotic behavior of the Green function (see the 
formula (l29l) V The mentioned factor produces that in the asymptotic limit G{r, r'; t) and dtG{r, r']t) decay in 
the form h{xi) where h{xi) is a sine or cosine function of xi- The behavior (|55)) is more complicated 

than the behavior (1361) or dazi), but usually we consider the initial conditions F{r',0) = 0 or dtF{r',0) = 0, 
that is, the behavior that we find in (l36l) or (|37|) . 

From our results we deduce that in the D-dimensional Minkowski spacetime the Dirac field decays with 
an inverse power of the time and this power is a half-integer for D-even and integral number for D odd. 
Furthermore we notice that for a fixed value of n, when the dimension of the spacetime increases, the Dirac 
field decays faster, since the power L + 3/2 of the time increases as D increases (see the expressions ([5])). 
In contrast to massless fields [23], [32], for the massive Dirac field our results point out that in Minkowski 
spacetime its late time behavior is not strongly affected whether D is even or odd, since in both cases we 
obtain a late time behavior of the form (IMl) or (|37]l . 

For the massive Dirac field propagating in four dimensional curved spacetimes, the previous results |15]- 
PP] show that its asymptotic behavior has two parts (see also [1]-[I1]). At the intermediate late time 
behavior the power law decay rate depends on the angular eigenvalue and the mass of the Dirac field, 
whereas in the very late time behavior the power law decay rate is of the form that is, it is independent 

of the parameters for the field and the spacetime m-m- In this work we find that in the D-dimensional 
Minkowski spacetime the asymptotic behavior of the massive Dirac field is different from the two previously 
described in four dimensional curved spherically symmetric spacetimes, since as t —>■ oo the power law decay 
rate of the massive Dirac field depends on the dimension of the spacetime, on the mode number n of the 
angular eigenvalues, and it does not depend on the mass of the Dirac field. Thus in curved spacetimes we can 
infer that for the massive Dirac field the dependence on the mass of the power law decay rate at intermediate 
late times is an effect of the curvature of the spacetime. 
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